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1.  INTRODUCTION 


Adiabatic  shear  bands  form  in  a  deforming  material  when  thermal  softening,  due  to  plastic 
heating,  is  stronger  than  work  hardening  and  rate  hardening  combined,  so  that  a  regime  with 
net  strain  softening  occurs.  Numerical  studies  by  Walter  (1992)  have  shown  that,  after  a 
period  of  relatively  slow  growth  of  perturbations,  it  often  happens  that  an  extreme  localization 
of  the  deformation  occurs  suddenly  so  that  subsequent  deformation  is  confined  almost  entirely 
to  a  narrow  strip  or  band  of  material,  which  becomes  very  hot  relative  to  the  adjacent,  almost 
nondeforming  material.  Within  these  narrow  bands,  often  measured  in  a  few  tens  of  microns 
or  less,  heat  production  from  plastic  working  is  nearly  balanced  by  heat  conduction  down 
extremely  steep  temperature  gradients. 

In  studies  on  adiabatic  shear  bands  (the  name  is  an  historical  misnomer),  a  recurring 
question  is  the  "width"  of  the  band,  which  should  be  definitively  furnished  by  the  theory.  It  is 
the  purpose  of  this  report  to  answer  that  question,  at  least  within  the  context  of  a  simple 
model,  by  examining  the  morphology  of  a  fully  formed  band  in  the  late  stages  of  deformation. 

2.  THE  MODEL 

Wright  and  Walter  (1987)  used  the  following  dimensional  equations  as  a  prototypical 
model  for  adiabatic  shear  banding: 


Sy  =  PV, 

(1.1) 

s,  =  n(vy  -  Yp) 

(1.2) 

pcA  =  k^  +  syp 

(1.3) 

s  =  (I  -  ad)  (1  +  bYp)m  . 

(1.4) 

Boundaries  at  y  =  ±h  are  assumed  to  be  moving  at  constant  velocity  and  to  be  fully  insulated. 
The  spatial  coordinate  perpendicular  to  the  slab  is  y,  time  is  t,  and  subscripts  in  y  or  t  indicate 
partial  derivatives.  Equation  1.1  represents  balance  of  linear  momentum  in  a  one-dimensional 
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slab  of  incompressible  material  undergoing  simple  shearing,  where  s  is  shear  stress,  v  is 
particle  velocity  in  the  plane  of  the  slab,  and  p  is  density.  Equation  1 .2  gives  the 
conventional,  additive  decomposition  of  elastic  and  plastic  strain  rates,  where  p  is  the  elastic 
shear  modulus,  vy  is  the  total  shear  strain  rate,  and  yp  is  the  plastic  strain  rate.  Equation  1.3 
represents  balance  of  energy,  where  d  is  temperature  (measured  from  an  arbitrary  reference 
level),  and  c  and  k  are  physical  constants  for  heat  capacity  and  thermal  conductivity, 
respectively.  It  has  been  tacitly  assumed  that  thermal  and  mechanical  influences  are  totally 
uncoupled  in  the  internal  energy  of  the  material,  that  the  dependence  of  internal  energy  on 
temperature  is  only  linear,  and  furthermore,  that  all  the  plastic  work  is  converted 
instantaneously  into  heat.  Equation  1.4  is  the  assumed  flow  law  for  an  elastic,  perfectly 
plastic  material  with  thermal  softening  and  rate  hardening,  where  Kq,  a,  b,  and  m  are  physical 
constants.  In  metals,  the  strain  rate  sensitivity  is  very  small,  m  <  .03  in  most  cases.  In  actual 
practice,  the  flow  law  is  used  as  a  definition  for  yp  when  s  >  Kq  (1  -  afl);  otherwise,  yp  =  0.  In 
this  simple  model,  work  hardening  has  been  ignored;  in  a  more  realistic  model,  an 
evolutionary  law  would  have  to  be  specified  for  the  strength  i^. 

A  numerical  study  by  Walter  (1992),  which  as  far  as  the  authors  are  aware,  contains  the 
only  extensive  treatment  of  fully  formed  shear  bands,  has  shown  that  the  model  may  be 
further  simplified  without  significant  loss  of  realism  in  certain  circumstances.  Thus,  if  the 
density  times  the  square  of  the  imposed  boundary  velocity  and  divided  by  a  characteristic 
stress  is  small  compared  to  1  (pvo/S0  «  1),  then  the  stress  stays  essentially  independent  of  y 
even  though  it  may  change  extremely  rapidly  with  time.  If  the  shear  modulus  is  large 
compared  to  a  characteristic  stress  (p/S0  »  1 ),  as  is  generally  true  for  metals,  then  after  some 
initial  transients,  the  plastic  strain  rate  and  the  total  strain  rate  become  nearly  identical. 

Finally,  since  byp  »  1  in  most  of  the  deformation,  the  flow  law  may  be  simplified  accordingly. 
When  these  simplifications  are  taken  into  account,  and  the  following  scheme  of 
nondimensionalization  has  been  used 

y  =  y/h.  T  =  it 

5  =  s/S0,  =  pcd/S0,  v  =  v/y0h 

P  =  P  (y,h)2/S0,  k  =  k/pcy,h2,  I  =  aSo/pc  (2) 
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then  a  reduced  set  of  equations,  now  in  nondimensional  form,  but  with  the  bars  dropped  for 
clarity,  is  as  follows: 


J* 

II 

o 

(3.1) 

a,  =  kd„  +  svy 

(3.2) 

s  =  (1  -aa)vy 

(3.3) 

The  nominal  applied  strain  rate  is  so  the  velocity  maintained  on  the  boundary  is  v0  =  -y0h, 
and  S0  as  iq,  (by0)m  is  a  characteristic  stress.  Boundary  conditions  for  the  reduced  problem  are 
taken  to  be  v(±1 ,  t)  =  ±1  and  dy(±1 ,  t)  =  0.  For  convenience,  only  symmetrical  solutions  will 
be  considered  further  (i.e.,  v(0,  t)  ■0,v(1,t)*1,  ay(0,  t)  =  ay(1,  t)  =  0). 

3.  A  NONLINEAR  EIGENVALUE  PROBLEM 

In  their  1987  paper,  Wright  and  Walter  noticed  that  in  their  numerical  solutions  the  strain 
rate  became  essentially  independent  of  time  in  the  late  stages  of  deformation,  although  the 
temperature  and  stress  continued  to  evolve.  Accordingly,  they  proposed  a  set  of  equations  to 
describe  the  terminal  behavior  (which  they  termed  quasi-steady)  but  they  did  not  analyse 
those  equations  further.  An  equivalent  set  of  equations  follows. 

First,  define  a  new  dependent  variable,  u(y),  from  vy  =  Au'1/m,  where  A  is  the  constant 
strain  rate  in  the  center  of  the  band,  and  u  does  not  depend  on  time.  From  the  flow  law 
(Equation  3.3),  u  =  Am(1  -  afl)/s,  and  from  Equation  3.1  and  Equation  3.2,  u  satisfies 

-  a  A1""  u^1*"”""  »  2  =  -a.  (4) 

u  s 

Since  it  has  been  assumed  that  u  depends  only  on  y,  and  s  only  on  t,  it  follows  by  the 
usual  separation  argument  that  a  is  a  constant.  Thus,  an  observation  about  numerical  results 
has  led  to  a  pair  of  ordinary  differential  equations. 
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uw  +  “u  -  1a1  *mu'1/m  -0,  6  +  as  -  0, 
k  k 


(5) 


the  solution  of  which  will  be  a  special  solution  of  the  partial  differential  equations  (Equations 
3.1 , 3.2,  and  3.3).  The  eigenvalue,  a,  and  the  normalizing  constant.  A,  must  be  determined 
as  part  of  the  solution. 

From  the  definition  of  u,  we  have  u(0)  =  1 ,  and  from  the  boundary  conditions  on  d,  we 
have  Uy(0)  =  uy  (1)  =  0.  Because  of  the  boundary  conditions  on  v,  the  integral  of  the  strain 
rate, 


1, 


implies 

j^u-1/mdy  -A"1. 


Then  integration  of  Equation  5,  together  with  the  use  of  this  result  and  the  boundary 
conditions  on  u,  implies 


i: 


udy  =  aA  m/a . 


(6) 


A  first  integral  of  Equation  5  that  satisfies  boundary  conditions  at  y  =  0  is 


u*  =  X(1  -  u-p)  -  v(u2  -  1) 


(7) 


where  the  constants  X,  p,  and  v  are  given  by 


i  _  2am 
k(1  -  m) 


(8.1) 
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p 


1  -  m 
m 


(8.2) 


v 


(8.3) 


Since  p  is  large  (p  >  30  for  metals),  the  phase  portrait  corresponding  to  Equation  7  looks 
nearly  circular  in  coordinates  uy  vs.  ^vu~  until  u  gets  close  to  1 ,  and  then  it  drops  sharply  back 
to  0,  as  shown  in  Figure  1.  The  radius  of  the  circle  is  nearly  A  +  v  ,  as  shown.  The  curve  is 
parameterized  by  the  y-coordinate  with  the  boundaries  y  -  0  and  y  =  1  lying  at  the 
intersections  where  \f  *  0.  Although  the  interval  0  £  y  £  1  may  correspond  to  multiple  circuits 
around  the  phase  diagram,  from  now  on  attention  will  be  focused  only  on  a  half  circuit  with 
u(0)  -  1  and  u(1)  -  A  +  X/v  ,  which  corresponds  to  the  sketches  for  u  and  vy,  as  shown  in 

Figure  1 .  Sketches  of  solutions  corresponding  to  multiple  circuits  may  be  obtained  by  obvious 
reflections  and  periodic  extensions.  It  should  be  noted  as  well  that  since  ad  =  1  -  sA^u,  the 
temperature  distribution  will  look  like  u,  only  linearly  scaled  and  turned  upside  down.  The 
important  point  to  note  is  that  the  temperature  distribution  will  appear  to  be  much  broader  than 
the  strain  rate  distribution,  which  will  be  intensely  concentrated  near  the  origin  due  to  the 
smallness  of  the  parameter  m. 

Figure  1  suggests  that  inner  and  outer  approximations  to  the  solution  should  be  sought 
corresponding  to  the  left  and  right  sides  of  the  phase  portrait.  In  the  outer  solution,  u  >  1  and 
u*«  1,  and  so  we  write 


(u'o)2  -  X  +  v  -  vuo  (9) 

with  boundary  condition  u'0(1)  =  0.  In  Equation  9,  at  this  stage  the  relative  size  of  X  and  v  is 
unknown.  The  solution  that  matches  the  boundary  condition  at  y  -  1  is 


cos>/v  (1  -  y). 


(10) 
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VvM- 


We  note  that  u0  -*  J  1  +  —  cos^v  +  JX  +  \  y  sinv/v  +  0(y2)  as  y  -»  0. 

For  the  inner  solution,  u  is  near  to  1 .  We  write  u(  =  1  +  q,  and  since  u  ?  -  e  for  q  «  1 , 
q  should  satisfy 

(q')z  =  X(1  -  e‘M)  -  v(2q  +  q2)  (11) 

As  q  approaches  0,  the  right-hand  side  of  Equation  1 1  approaches  (Xp  -  2v)q.  Since  p  is 
large  compared  to  1,  even  if  v  and  X,  are  comparable  in  size,  it  still  turns  out  that  2v  «  Xp,  so 
the  second  term  on  the  right-hand  side  of  Equation  1 1  may  be  dropped  for  now,  subject  to 
verification  of  orders  of  magnitude  later.  Putting  r  =  e'™,  by  quadrature  we  find  the  solution  to 

(a  —  \ 

Equation  11  to  be  r  =  sech2  — pyX,y  ,  or  in  the  original  variables, 

l2 

of  1 

Ui  *  1  +  _ln  cosh— 

P  l  2 

We  note  that  u.  -»  1  -  —\n2  +  >/Xy  as  y  -»  » . 

P 

To  close  the  phase  portraits,  we  will  try  to  patch  the  two  solutions  together  at  a  common 
point.  That  is,  find  a  point  y  =  ym  such  that 

ut(y J  =  MyJ,  u/  (ym)  =  uo  (ym)  (13) 

or,  in  detail, 

cosy^O  -  ym)  -  1  +  -  In  cosh  Ip  JXym  (14.1) 

P  \  2 

rT7V  sin/7 (1  -  yj  =  v^tanh^p^y,,, 


(14.2) 


where  higher  order  terms  in  ym  and  p*1  have  been  neglected.  Within  terms  of  the  order 
(pVX)-1  ,  the  linear  terms  in  ym  are  eliminated  by  choosing 


s\njv  «  .  v^~  -  . 
A  +  v 


(18) 


Equation  18  could  have  been  obtained  directly  by  matching  the  slopes  of  inner  and  outer 
solutions,  i.e.,  u',  (°°)  *  u'0  ( 0 ) . 

It  is  clear  that  there  are  an  infinity  of  solutions  for  v  as  a  function  of  X  since  the  right-hand 
side  of  Equation  1 8  decreases  monotonically  from  1  towards  0  as  v  increases,  whereas  the 
left-hand  side  is  periodic  with  repeated  maxima  of  1 .  Thus,  there  will  be  two  solutions  of 
Equation  18  for  each  maximum  in  the  sine  function.  Since  X  is  expected  to  be  large,  the  right- 
hand  side  decreases  slowly  with  v,  and  the  first  intersection,  which  is  also  consistent  with  the 
desired  phase  portrait,  will  occur  near  the  first  maximum.  It  turns  out  that 

=  £  -  e  where  e  *  0(X~1/2)  (19) 

and,  consequently,  the  decay  rate  of  the  stress  is  given  by 

a  =  flk.  (20) 


Since  Equation  19  shows  that  v  =  0(1),  it  would  have  been  justified  to  ignore  v  in  comparison 
with  X  in  Equation  9. 

In  principle,  from  the  inverses  of  Equations  16  and  15.2,  we  could  now  go  back  to  find  ym, 
which  from  Equation  15.2  must  be  0{pAV’  •  But  rather  than  doing  that,  we  take  the 
thickness  of  the  boundary  layer  to  be 
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-2/p/T-f. _ SS _ Y- 

^(1  -  mJaA1  *m > 

where  the  argument  of  cosh  in  Equation  12  is  1. 

The  amplitude  A  and  band  width  8  may  now  be  determined.  The  left-hand  side  of 
Equation  6  may  be  evaluated  approximately  as  follows: 

J^udy  -  J\,dy  .  £  u„dy 

-  8  +  v-1/2;i  +  (X/v)  [sinv/v  (1  -  8)] 

-  v-^(1  +0 (X*,/2)).  (21) 

Note  that  the  value  of  the  integral  comes  all  from  the  outer  solution  for  8  «  1 .  Since  the  right 
hand  side  of  Equation  6  is  aAm/kv  after  using  Equation  8.3,  we  have  fk  =*  aAm/k,  or  with  the 
aid  of  Equation  8.1, 


fl-m  a 

<  2  mkj 


(22) 


It  is  now  easy  to  work  out  that 


5  =  1/A. 


(23) 


Clearly,  since  m  «  1 ,  the  condition  that  8  «  1  and  A  »  1  reduces  approximately  to 
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Also,  since 


2 v 
Xp 


* 

2m  }■ 

4 

[1  -  mj 

1  ♦  m 


the  condition  that  was  required  to  justify  the  inner  solution,  2vAp  «  1 ,  reduces  approximately 
to  k/a  «  m'1/2,  which  is  more  stringent  than  the  previous  condition.  However,  the  most 
stringent  condition  of  ail  comes  from  the  condition 


»  1, 


which  was  required  to  obtain  Equation  19.  Since  the  first  factor  is  0(1)  for  small  m,  the 
approximate  requirement  is  that  (k/a)2  «  1 .  From  the  nondimensionalization  in  Equation  2,  it 
may  be  seen  that  the  strong  inequalities  fail  if  the  nominal  applied  strain  rate,  is  too  small. 
Finally,  by  back  substitution  we  see  that 


x(lnx)1* 


n'T-rm 

'1 

-  m> 

i  *  m 

£ 

Z<1  -  m) 

a 

n 

m  J 

1 

-  m 


which  tends  to  infinity  as  m  tends  to  zero  (p  tends  to  infinity),  and  thus  justifies  the 
assumptions  associated  with  Equations  15, 16,  and  17. 

When  Equations  22  and  23  are  expressed  in  dimensional  terms,  it  becomes  clear  that  the 
maximum  strain  rate  in  the  shear  band  and  the  thickness  of  the  shear  band  depend  only  on 
the  imposed  velocity  difference  across  the  band,  rather  than  on  the  imposed  strain  rate  and 
slab  width  separately.  Thus,  converting  Equation  22  back  to  dimensional  quantities,  we  have 
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Tm«  *  To^  “  7o 


(  \1/(1  -"»> 

1  -  m  aS0  pcv0h2 


2m  pc  k 

f  v/(i-m) 


) 


1  -  m  K  m 


2m 


In  the  final  form  of  Equation  24,  the  expression 


2/(1  -m) 
v0 


(24) 


a  =  ^1  iH  =  —  (by. )“ 

S0  3d  S0  3d 


has  been  used,  and  3o0/3d  represents  the  rate  of  thermal  softening  under  quasi-static 
conditions  and  ambient  temperature.  In  the  present  model,  a0  =  k^,  (1  -  ad),  and  the 
derivative  is  evaluated  at  d  =  0.  Similarly,  the  dimensional  form  of  Equation  23  becomes 
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h/A 


( 

1  -  m  -9o0/3d 
2m  k 


1  ♦  m 


j 


(25) 


4.  NUMERICAL  EXAMPLE  AND  CONCLUSIONS 

Wright  and  Walter  (1987)  used  the  following  nondimensional  data  in  their  full  finite  element 
calculations  with  x  =  500s'1: 

m  =  0.0251 ,  k  =  0.0022,  a  =  0.1 04  x  (5  x  1 06)oo2S1. 

The  softening  coefficient,  a,  has  been  multiplied  here  by  the  factor  (bY0)m  to  account  for  the 
fact  that  Wright  and  Walter  scaled  stresses  by  Kq,  whereas  stresses  have  been  scaled  by  S0  in 
this  report.  The  data  lead  to  the  values  A  =  1 ,628  and  log10  5  =  3.212,  which  appear  to  be 
extremely  accurate  when  compared  to  Wright  and  Walter’s  Figure  7  (see  Figure  2  of  this 
report).  Note  that  the  two  curves  differ  only  slightly  and  only  at  the  top;  the  estimate  for  5  lies 
at  the  point  where  the  strain  rate  has  fallen  to  approximately  40%  of  its  maximum  value. 
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Figure  2.  Comparison  of  Present  Asymptotic  Solution  With  Finite  Element  Solution  Due  to 
Wright  and  Walter  (19871. 
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Equations  24  and  25  give  the  basic  scaling  laws  for  the  intensity  and  width  of  the  strain 
rate  profile.  For  m  «  1,  as  in  metals,  neither  ^  nor  the  dimensional  value  of  8  is  sensitive  to 
the  exact  value  of  b.  As  expected,  the  maximum  strain  rate  increases  with  the  rate  of 
softening,  but  decreases  with  strain  rate  sensitivity  and  thermal  conductivity,  and  the  width  of 
the  strain  rate  profile  varies  in  the  opposite  sense.  As  has  already  been  noted,  both  quantities 
depend  only  on  v0,  and  not  on  ^  and  h  separately.  It  is  interesting  to  note  that  neither 
quantity  depends  on  either  the  density  or  the  heat  capacity,  which  is  a  consequence  of  the 
nearly  steady  feature  of  fully  formed  bands. 

One  last  consequence  of  Equations  24  and  25,  which  may  have  some  experimental 
significance  even  though  it  has  been  obtained  with  a  special  flow  law,  is  that  the  product  of 
the  maximum  strain  rate  and  the  band  width  is  equal  to  the  imposed  velocity 

-  Vo  •  (26) 

The  width  5  as  obtained  in  this  analysis  is  the  distance  from  the  center  of  the  band  at 
which  the  strain  rate  is  less  than  the  maximum  value  by  a  factor  of 


(or  about  2.5),  but  from  there  on  the  strain  rate  decreases  so  rapidly  that  this  estimate  of  width 
should  be  as  good  as  any. 

In  an  experimental  sense  the  "width"  of  a  shear  band  is  often  associated  with  the  width  of 
a  "white  etching  band,"  which  in  the  older  literature  has  often  been  referred  to  as  a 
"transformed  band,"  implying  that  a  metallurgical  transformation  has  occurred  that  causes  the 
material  in  the  band  to  etch  in  a  manner  different  from  the  surrounding  material.  If  that 
"transformation"  is  associated  only  with  the  extreme  straining  in  the  fully  formed  band,  leading 
perhaps  to  ultra  fine  grain  refinement  with  little  or  no  recrystalization,  as  recent  work  suggests 
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for  4340  steels  (e.g.,  Beatty  et  al.  1990*),  then  estimates  of  the  kind  given  in  Equations  25 
and  26  should  be  accurate  and  useful  for  "transformed"  bands,  as  well  as  for  the  simpler 
"deformed"  bands,  where  localization  has  occurred  without  a  change  in  etching  behavior.  On 
the  other  hand,  metallurgical  transformations  have  not  yet  been  included  in  any  theoretical 
treatment  and,  in  any  case,  may  be  associated  with  the  thermal  profile,  which  is  much  wider 
than  the  strain  rate  profile.  Therefore,  a  more  comprehensive  estimate  of  "band  width"  for 
cases  with  true  transformations  remains  lacking. 

As  a  final  comment,  the  reader  is  reminded  that  the  results  in  this  paper  are  specific  for 
the  assumed  flow  law  in  Equation  3.3,  which  is  intended  to  describe  a  perfectly  plastic 
material  with  linear  thermal  softening  and  power  law  rate  hardening.  In  fact,  the  separation  of 
variables  that  is  shown  in  Equation  5  only  occurs  exactly  in  the  case  of  linear  softening. 
Nevertheless,  the  methods  and  results  shown  here  are  broadly  suggestive  for  more  general 
cases,  for  example,  when  the  softening  rate  varies  slowly  with  temperature,  or  when  work 
hardening  nears  saturation. 


*  The  authors  are  grateful  to  a  reviewer  for  calling  this  reference  to  their  attention. 
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